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ABSTRACT 

We present high-dynamic-range images of circumstellar dust around HR 4796A that were obtained 
with MIRLIN at the Keck II telescope at A = 7.9, 10.3, 12.5 and 24.5 /zm. We also present a new 
continuum measurement at 350 /im obtained at the Caltech Submillimeter Observatory. Emission is 
resolved in Keck images at 12.5 and 24.5 ^m with PSF FWHM's of 0.37" and 0.55", respectively, and 
confirms the presence of an outer ring centered at 70 AU. Unresolved excess infrared emission is also 
detected at the stellar position and must originate well within 13 AU of the star. A model of dust 
emission fit to flux densities at 12.5, 20.8, and 24.5 /im indicates dust grains are located 4^ AU from 
the star with effective size, 28±6 /im, and an associated temperature of 260±40 K. 

We simulate all extant data with a simple model of exozodiacal dust and an outer exo-Kuiper ring. 
A two-component outer ring is necessary to fit both Keck thermal infrared and HST scattered-light 
images. Bayesian parameter estimates yield a total cross-sectional area of 0.055 AU 2 for grains roughly 
4 AU from the star and an outer-dust disk composed of a narrow large-grain ring embedded within 
a wider ring of smaller grains. The narrow ring is 14±1 AU wide with inner radius 66±1 AU and 
total cross-sectional area 245 AU 2 . The outer ring is 80±15 AU wide with inner radius 45±5 AU and 
total cross-sectional area 90 AU 2 . Dust grains in the narrow ring are about 10 times larger and have 
lower albedos than those in the wider ring. These properties are consistent with a picture in which 
radiation pressure dominates the dispersal of an exo-Kuiper belt. 

Subject headings: circumstellar matter - infrared: stars - planetary systems: formation - planetary sys- 
tems: protoplanetary disks - solar system: formation - star: individual(HR 4796A) 



1. INTRODUCTION 

In the last half decade, advances in high-resolution 
long-wavelength techniques have yielded an increased 
number of spatially resolved images of thermal emis- 
sion from dusty disks around nearby stars (Greaves et 
al. 1998; Holland et al. 1998; Jayawardhana et al. 1998 
[Jay98]; Koerner et al. 1998 [K98]; Holland et al. 1999; 
Koerner, Ostroff & Sargent 2001; Wilner et al. 2002; 
Wahhaj et al. 2003; Weinberger et al. 2003). These 
typically reveal dust masses located at several tens of 
AUs from the central star in a region analagous to the 
Edgeworth-Kuiper Belt in our own solar system. In most 
instances, emission is markedly reduced close to the star. 
This inner clearing was previously inferred in some cases 
on the basis of the spectral energy distribution (SED) 
(e.g., Jura 1998). Dust located here is most easily de- 
tected at mid-infrared wavelengths where stellar photo- 
spheric radiation is relatively diminished and circumstel- 
lar dust emission is at a maximum. If located within an 
ice-condensation radius, it can be considered to be anal- 
ogous to the solar system's zodiacal dust, a population 
of grains that emanates from asteroid collisions and the 
disintegration of comet tails (Dermott et al. 1994). The 
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presence of such an "exozodiacal" dust component was 
inferred by K98 on the basis of high-resolution (~ 0.4") 
12 and 20 jum images of HR 4796A. 

The detection of exozodiacal dust is of considerable 
importance to efforts aimed at the detection of earth- 
like planets (Backman et al. 1998). In light of a rapid 
dispersal timescale (< 10 4 yr; Gustafson 1994), current 
detections imply an origin in the collision of minor plan- 
etary bodies and perhaps resonant trapping by larger 
ones. The dust also may impact planetary detection 
techniques as a source of confusing radiation (Backman 
et al. 1998). Despite their importance, detections are 
currently few in number. Evidence for exozodiacal dust 
around HR 4796A was recently brought into question by 
Telesco et al. 2000 (TOO) and Li & Lunine 2003 (LL03). 
Resolution of the issue is complicated by the fact that 
a prominent outer dust ring provides a confusing signa- 
ture at the resolution of HR 4796 A images. Here we 
present additional high-resolution images of HR 4796A 
with modeling designed to better understand the radial 
structure of its circumstellar dust, including whether or 
not an exozodiacal dust component is required by the 
data. 

HR 4796A is an A0V star (Houk & Sowell 1985) lo- 
cated at a distance of 67±3 AU (Hipparcos catalog) with 
an age of 8±2 million years (Stauffer et al. 1995). K98 
and Jay98 first resolved the disk in thermal infrared im- 
ages. Jay98 noted that a line-cut through their 18 fim 
image could be simulated with a ring 50 to 110 AU in 
extent. K98 simulated Keck images at higher resolu- 
tion and demonstrated over a large range of models that 
an inner clearing was demanded by the data. The most 
probable values of disk properties yielded inner and outer 
radii of 55 and 80 AU, respectively, for a disk inclined 
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~ 73° from face on. K98 also found excesses at the stel- 
lar position at 12 /im and 20 /mi with values that implied 
a dust temperature of ~ 250K and, in turn, a location 
~ 4.5 AU from the star. Schneider et al. 1999 (S99) im- 
aged scattered light from the HR 4796A outer ring with 
HST/NICMOS and confirmed the presence of a narrow 
(17 AU) outer ring centered at 70 AU. TOO and Wyatt 
et al. 1999 (W99) analyzed Keck images at 10 /urn and 
18 /im and derived a radial surface density to the outer 
ring that rose exponentially from ~ 45 AU, doubled be- 
tween 60 and 70 AU, and fell like r" 3 outside 70 AU with 
an outer cutoff at ~ 130 AU. Grain sizes were estimated 
to be 2-3 /im in contrast to 30 /im grains found by K98. 
TOO also found that any excess detected in their images at 
the stellar position was consistent with their photomet- 
ric uncertainties. Augereau et al. 1999 presented models 
of the SED of HR 4796A which predicted dust within 
10 AU of the star. LL03 showed that the SED could 
be modeled without a zodiacal dust component. Sitko 
et al. 2000 observed that the silicate emission feature is 
very weak, suggesting that there aren't many small sil- 
icate grains near the star. He also estimated that only 
half the flux in the 8-13.5 /im range is photospheric, de- 
tecting more excess than in the 12.5 /im MIRLIN image 
in K98. Here we present new data and model all extant 
measurements simultaneously. These allow us to resolve 
disputes about the properties of dust around HR 4796A. 

2. OBSERVATION 

HR 4796A was observed on UT 11 k 14 June 1998, 30 
January 1999, and 3 February 1999 with the mid-infrared 
camera MIRLIN at the Keck II telescope. MIRLIN's de- 
tector is a Boeing 128x128 pixel, high-flux Si: As BIB 
array with a Keck II plate scale of 0.138" per pixel and 
17.5" field of view (Ressler et al. 1994). Chopping with 
the Keck II secondary mirror was carried out in a North- 
South direction at a 4 Hz rate with a throw of 7", In 
between chop sequences, the telescope was nodded East- 
West with a similar throw. Since thermal IR emission 
from HR 4796A is confined to within a maximum length 
of 5", the 7" chop/nod throw was sufficient for differenc- 
ing of any pair of images without overlap. Observations 
were carried out in filters centered at A = 7.9, 10.3, 12.5 
and 24.5 /im with widths 0.76, 1.01, 1.16 and 0.76 /im and 
for on-source integration times of 3, 6, 6, and 52 minutes 
respectively. Images on doublc-diffcrcnccd frames were 
cross-correlated with a point spread function (PSF) to 
determine a first estimate of the centroid position and 
were then coadded. The resulting image was then used 
in place of the PSF to repeat the process with improved 
centroiding and produce a new stacked image. This pro- 
cedure was repeated multiple times until there was no 
perceptible change in the final coadded image on suc- 
cessive iterations. The final image was sized to 64 x 
64 pixels (8.8" x 8.8") to optimally window the emission 
from HR 4796A and adjacent background. Infrared stan- 
dards (3 Leo, 2 Cen, a Sco and a Boo were observed in 
the same way at 7.9 (and 10.3), 12.5, 20.8 and 24.5 /im 
respectively at similar airmasses. 

To better constrain the properties of grains around 
HR 4796A, a new continuum measurement was obtained 
on UT May 23, 1998 at 350 /im with the SHARC lin- 
ear bolometer array camera at the Caltech Submillimetcr 
Observatory. Atmospheric conditions were good with a 



steady T225 of 0.36. The flux calibrator was Uranus with 
an angular size of 3.59" and a flux density of 242 Jy 
on that date. The resulting 350 /im flux density of 
HR 4796A was observed to be 160±42 mJy and is listed 
in Table. 1. 

Images of HR 4796A are unresolved at A = 7.9 and 
10.3 /im with respective PSF FWHM of 0.48" and 
0.45" and are not displayed here. Flux densities were 
calculated by comparison with images of (3 Leo ob- 
tained shortly after those of HR 4796A. Flux densities 
of 307±44 mJy and 218±24 mJy at 7.9 and 10.3 /im 
were derived, assuming values of 10.68 Jy and 6.44 Jy 
for (3 Leo at the same wavelengths (Koerner et al. 2000). 
Photospheric estimates 283 mJy at 7.9 /im and 168 mJy 
at 10.3/im, were derived by the fit of a model A0V star 
(Kurucz 1993) to measurements of HR 4796A at wave- 
lengths shorter than 3 /im. It seems that there is a small 
excess at 10.3 /im and an even smaller excess or possibly 
none at 7.9 /im. 

The 12.5 /im image is shown in Fig. [Ha). The 64 x 64 
image was finely gridded to 256 x 256 pixels and then 
smoothed by convolution with a circular hat function of 
radius equal to the FWHM of a Boo (PSF star) which 
was observed immediately after HR 4796A. The resolu- 
tion of the image is 0.37". Peak brightnesses are calcu- 
lated as the flux in the brightest pixel, divided by the 
pixel area. The central peak brightness was measured 
to be 723±44 mJy/arcsec 2 while that expected from the 
photosphere is 387 mJy/arcsec 2 . The total flux at 12.5 
fim is 231±14 mJy, while the photospheric flux at 12.5 
/im is 114 mJy. The long axis enclosed by the 3<r contour 
is 2.1" (141 AU). NE and SW peaks are detected with 
total separation 1.9" (127 AU). The NE and SW peaks 
are not equal in brightness, but have a flux ratio of 1.13. 

The 20.8 /im image displayed in Fig. db) is a re- 
reduced version of data presented in K98. Shifting and 
adding was carried out as described above, starting with 
a model of the emission derived from earlier versions of 
the image. The result exhibits improved signal to noise 
and a flux density, 1.62±0.16 Jy, that is consistent with 
the previously reported 1.88±0.17 Jy. The final shift- 
and-add image was smoothed to match the 0.4" resolu- 
tion of the other images in Fig.^ Emission peaks NE and 
SW of the star are separated by 1.85" (124 AU) in the 
final smoothed image. The central peak is brightest and 
0.93" (62 AU) from the SW peak and the same distance 
from the NE peak. Peak brightness values are 629, 668, 
and 623 mJy/arcsec 2 at the NE, center, and SW posi- 
tions, respectively. The photospheric flux at A= 20.8 /im 
is 41 mJy, and which corresponds to a peak brightness 
of atmost 82 mJy/arcsec 2 . 

The 24.5 /im image is shown in Fig. ^c) and was re- 
duced in the same way as the 12.5 /im image with a Boo 
used as the psf star. Image resolution was 0.55". The 
expected flux from the photosphere is 30 mJy, while the 
total measured flux for HR 4796A was 2100±170 mJy. 
The expected photospheric brightness in the peak pixel 
is at most 44 mJy/arcsec 2 . Comparing this to the cen- 
tral peak brightness of 585 mJy/arcsec 2 it is immediately 
obvious that most of the emission in the central pixel is 
non-photosphcric and probably the result of hot circum- 
stellar dust. The emission extends 1.68" (113 AU) from 
the NE peak to SW peak. The extent of the 3tr contour 
is 3.6" (242 AU). The central peak is 0.74" (50 AU) from 
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Fig. 1. — KECK/MIRLIN contour maps of HR 4796A from new observations at thermal infrared wavelengths and the 
HST/NICMOS image from S99. North is aligned with the up position. Contour images of the PSF stars, 2 Cen, a Sco and 
q Boo, are displayed in panels a through c as insets with contour intervals of 10%, starting at the 60% level. The panel d-inset 
shows a TINYTIM PSF at the resolution (0.12") of the HST/NICMOS observation at A = l.l^m. The contour levels for the 
d-inset are the same as in the other panels, (a) Emission at A = 12.5 /im . Lowest contour is at the 2a level (32 mjy/arcsec 2 ). 
Higher contour levels are at 2a x (io 068n ) for the n th contour, (b) Emission at A = 20.8 /im with la (47 mjy/arcsec 2 ) contour 
intervals and contours starting at the 2a level, (c) Emission at A = 24.5 /jm contoured as in b) with la = 39 mjy/arcsec 2 . (d) 
Scattered light at A = 1.1 /jm . This observation of HR 4796A was made with HST/NICMOS (S99) at 1.1/an . The NICMOS 
coronagraph was used to mask the star at the center of the image. Scattered light within 0.6" of the star was deemed artifactual 
and is not shown in this image. The lowest contour is at 2a ( 1.8 mjy/arcsec 2 ) and successive contours are drawn at 2a intervals. 



the SW peak and 0.94" (63 AU) from the NE peak. Peak 
brightnesses are 637, 617, and 585 mjy/arcsec 2 at the 
NE, SW and central positions, respectively. The NE/SW 
brightness ratio is 637:617 = 1.04. 

A very high-resolution image of light scattered by the 
HR 4796A ring was obtained by S99 with the NICMOS 
coronagraphic camera on HST. We have contoured the 
1.1 £im image and displayed it in Fig.^d) for comparison 
to the thermal infrared images. The region within 0.6" 
of the star has artifactual emission from coronagraphic 
scattering and had been masked off. At a resolution of 
0.12" (8 AU), the emission is extended 2.06" (138 AU) 
from NE to SW peak, and exhibits a brightness asymme- 
try (NE:SW) of 14.4 : 12.2 = 1.18. The extent of the 2a 
contour line is much more confined than in the thermal 



infrared images. 

The images in Fig. ^confirm the general properties of 
ring emission inferred from previous high-resolution im- 
ages (K98, S99, TOO). Flux density measurements at the 
stellar position also appear to substantiate earlier evi- 
dence for a hot dust component close to the star (K98). 
The latter interpretation requires careful consideration, 
however, since the highly inclined ring passes within 0.3" 
of the star, and the PSF FWHM is 0.4". It is necessary 
to carefully account for the contribution of ring emis- 
sion at the stellar position to identify the presence of hot 
dust near the star with certainty. In addition, the struc- 
ture of the ring, itself, may be better ascertained with 
detailed analyses that combine information from all im- 
ages. To accomplish both goals, we compare a variety of 
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TABLE 1 
Flux Densities for HR 4796A 
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ring models to the Keck/MIRLIN and HST images. 

3. MODELING 

Our modeling approach is designed to best interpret 
all available observations of HR 4796A, including op- 
tical and thermal infrared images and flux densities at 
longer wavelengths. As a first step, we examine individ- 
ual Keck/MIRLIN images for evidence of dust close to 
the star. Then we quantitatively compare models of dust 
emission by fitting them simultaneously to all available 
data. To successfully deduce key disk properties in this 
manner, it is necessary to strike a balance between the 
need for a detailed physical description (including grain 
shape, size and composition) and an appropriate num- 
ber of free parameters. To best achieve this compromise, 
we use a simple model of emission from a flat optically 
thin disk as described in Backman, Gillett, & Witteborn 
(1992). It will be evident from the modeling process that 
further complications in disk structure, such as a flared 
disk, or a disk with a thick inner edge, etc., are not re- 
quired at the resolution of extant data. According to 
the Backman et al. model, thermal radiation from an 
annulus of width dr at distance r from the star is 



f t (r) = (I - u>)a(r)e x B[T p (r), A] [~^-srj JV, 

where <j(r) is the fractional surface density, T p (r) is 
the grain temperature, and D is the distance to HR 
4796A, 67 pc. The optical depth is defined as r(r, v) 
= (I — uj)a(r)e\, where u) is the albedo. For moder- 
ately absorbing dielectrics and an effective grain radius, 
a, the radiative efficiency is e\ — 1.5a/ X for A > 1. 5a 
and e\ = 1 for A < 1. 5a (Greenberg 1979). The grain 
temperature for efficient absorbers and inefficient emit- 
ters is T p (r) = 468( J L*/A ) - 2 (r/lAU)- " 4 as calculated 
from radiative balance equations (Backman & Paresce 
1993). is the stellar luminosity in solar units. The 



fragmentation process caused by inelastic collisions be- 
tween large grains results in a steady state size distri- 
bution described by n(a) ~ a~ 3 - 5 (Dohnanyi 1969). It 
has been shown that the spectral energy distribution re- 
sulting from such a collisional cascade of grains can be 
approximated by the SED of a population of grains of a 
single size (Backman, Gillett, & Witteborn 1992). The 
minimum grain radius in the distribution is a m i n ~ a/4, 
where a is the effective radius of the representive grains. 

To simulate the HST/NICMOS image from S99 in 
Fig.^d), we use the following prescription for an emis- 
sion model. The scattered-light emission from an annulus 
of width dr at distance r from the star is 

where .B* and R+ are the surface intensity and radius of 
the star. It is important to note that grain size does not 
affect the intensity of scattered-light emission, as grain 
temperature is not a factor. Therefore, the grain radius, 
a will not be a parameter in scattered-light modeling. 

In keeping with simple theoretical expectations, mod- 
els used to interpret low-resolution observations of debris 
disks often approximate the surface density, cr, as propor- 
tional to r -7 , where r is the stellocentric distance (Gillett 
1986; Diner & Appleby 1986; Nakano 1988; Artymowicz, 
Burrows & Paresce 1989 and Backman, Gillett & Wit- 
teborn 1992). High-resolution imaging has revealed a 
ring-like structure in many cases, however, necessitating 
the specification of other model parameters such as inner 
and outer radii Ri n and R ou t and ring width AR. 

To estimate photospheric emission from HR 4796A, we 
scale a model of the emission from an A0V star by Ku- 
rucz (1993). Model fitting uses all available photometric 
measurements of HR 4796A at wavelengths shortward of 
3 /im. To simulate images, the resulting photospheric 
flux density is used to scale a point source at the stellar 
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Fig. 2. — Models of stellar emission plus thermal emission from an outer ring scaled to match the outer disk's intensity are shown at 
12.5, 20.8 and 24.5 fim. The solid lines represent line cuts through data. The dashed lines represent line cuts through thermal-emission 
models of HR 4796A with a simple outer-ring as described in the text. All line cuts have a PA of 26° and pass through both the NE and 
SW peaks of emission. 



position. As a first approximation, we then add emission 
from a single-ring model and convolve the entire image 
with an appropriate PSF star. As demanded by the data, 
we then consider the addition of multiple rings in a sim- 
ilar modeling process. 

To calculate the range of acceptable values for model 
parameters as dictated by a simultaneous comparison to 
all data, we use a Bayesian approach that assigns a rel- 
ative probability, e -1 / 2S,,x ™, to each model. Here, \n 2 
refers to x 2 as defined for the nth data set. This requires 
a consistent way of weighting x 2 from the image fits with 
those from SED fits. Thus, for each image, all pixels of 
intensity above 3a are counted and divided by the area 
(in pixels) enclosed by the half-maximum contour of the 
corresponding PSF. This gives a measure of the number 
of data points represented by an image. The x 2 s from 
each of the image fits is normalized using this number. 

Probabilities are calculated for models varying over a 
large parameter space, and the relative probability of a 
specific parameter value is derived by summing proba- 
bilities for all models with that value. In this way, prob- 
ability distributions are built up for the values of each 
parameter (e.g. Lay, Carlstrom, & Hills 1997). As a 
prelude to modeling the entire data set, we first consider 
the question of whether or not the images mandate the 
presence of material close to the star that is analogous 
to the Solar System's zodiacal dust. 

3.1. Exozodiacal Dust 

To date, consensus is lacking in published interpreta- 
tions of evidence for an inner zodiacal dust component 
to the HR 4796A system. K98 derived an infrared ex- 
cess at the stellar position by fitting a scaled PSF to 
high-resolution Keck images at 12.5 and 20.8 /*m and 
comparing the results to expected photosphcric levels. 
TOO failed to detect this excess at 18.2 /im, however, by 
carrying out a similar analysis on Keck images taken un- 



der non-ideal photometric conditions. LL03 were able to 
fit the Spectral Energy Distribution (SED) without an 
inner dust component by varying grain composition and 
using outer-disk constraints supplied by the 1.1 /iin HST 
image (S99). 

The contradiction implied by these results reflects the 
difficulty of the analysis. Thermal infrared flux cali- 
bration is difficult from the ground, and the inner dust 
component is at the limit of available spatial resolution. 
For example, TOO derived a total flux density 20% lower 
than their previous measurement in the same filter at a 
different telescope (Jay98). If the latter measurement 
is in error, it would lead to an underestimate of the 
point-source excess after subtraction of flux from a pho- 
tosphcric model. It is also not surprising that LL03 did 
not find the case for inner dust compelling on the basis of 
SED shape alone, since the function that corresponds to 
emission at a given stellocentric distance is quite broad. 
The dust component inferred by K98 produces a negli- 
gible contribution to the overall SED and can be prop- 
erly substantiated only by high-resolution imaging data. 
Keck thermal infrared images of HR 4796A easily sep- 
arate emission at the stellar position from that of the 
ring along the long axis of emission (8 ~ 1"). However, 
the ring is oriented 27° from edge on, and the observa- 
tions have insufficient resolution across the minor axis 
(9 < 0.25"). Consequently, simple PSF fitting will be 
subject to an uncertainty from contamination by outer- 
ring emission. Here, we present a modeling approach 
that accounts for confusion from the outer ring in order 
to interpret new high-resolution images. 

To make an estimate of the contribution of outer-ring 
radiation to flux detected at the stellar position, we first 
fit a simple model of outer-ring emission solely to the re- 
gion outside of 35 AU (0.52") from the star. Our initial 
model was a single flat ring with density er(r) ~ r~ 7 , 
rotational axis inclined 73° to the line of sight (S99), 
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and with the long axis aligned at PA 26° . A point 
source with photospheric flux density was added at the 
stellar position and the entire model image convolved 
with a PSF star. Ring parameters were varied in inde- 
pendent fits to each of 3 images at A = f 2.5, 20.8, and 
24.5/im. Fig. 0shows line cuts through both the images 
and best-fit models. It is immediately evident that, in 
each case, these models fail to reproduce an adequate 
amount of emission at the stellar position. This compar- 
ison strongly supports the inference of an inner zodiacal 
dust component as concluded by K98. 
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Fig. 3. — Photometry of HR 4796A from Table 1 are denoted 
as asterisks. Zodiacal excesses calculated from point-source fitting 
to images at A = 12.5, 20.8 and 24.5 /jm are shown as circles at 
the bottom-left corner of the figure. Estimated flux densities from 
the zodiacal component at 10.3 and 7.9 fim are shown as "x". The 
dashed line represents the simulated SED of a model fitted to the 
zodiacal excesses. This model was a dust ring (T ~ 280 K) with 
radius 3.3 AU, width 0.5 AU, grain radius 14 fim and a = 0.0045. 



We make a quantitative estimate of the contribution 
of emission from inner dust by re-fitting our model to 
unmasked versions of the thermal infrared images while 
varying the model flux density at the stellar position. 
Expected photospheric contributions of 114, 41, and 30 
mJy at 12.5, 20.8, and 24.5 /tm , respectively, were 
subtracted from the best-fit result to yield estimates of 
92±41, 150±38 and 87±22 mJy for the inner component 
flux densities. The results are plotted in Fig. 03 It is 
evident from Fig. [3] that the inner component of dust 
contributes neglibly to the total SED at all wavelengths. 
Hence, the above result is consistent with the null result 
of LL03. As a first approximation to estimate physi- 
cal properties of a zodiacal dust component, we assume 
grains are located in a single ring of width 0.5 AU and 
vary the ring inner radius, r in , surface density, a, and ef- 
fective grain radius, a. Best-fit values are r in — 4^ AU, 
g = 0.0046±0.001, and a = 14±3 /im with uncertainties 
determined by the Bayesian approach described earlier. 
This implies an average grain temperature of 260±40 K. 
Ice grains of this size in a gas-free environment would 
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Fig. 4. — Probability distribution for the radius of the zodiacal 
ring obtained by the Bayesian approach. The best estimate for the 
radius is 4+j AU. 



have to be at a radial separation of atleast 43 AU to 
survive sublimation for any appreciable amount of time 
(Isobe 1970). To evaluate the consistency of our 10.3 and 
7.9 /im photometry with this simulation, we also plot the 
estimated contribution to the inner component at these 
wavelengths in Fig.[3J Assuming that the outer ring SED 
can be approximated by a 110K blackbody (Jura 1998) 
and that at 12.5 /im the outer ring contribution is ~ 25 
mJy, we extrapolate the outer ring contribution to be ~ 
5 mJy at 10.3 /im and and 0.2 mJy at 7.9 /im. Thus 
the inner component flux density is 45±24 mJy at 10.3 
/tm and 24±44 mJy at 7.9 /tm. From Fig. [31 it is evi- 
dent that our zodiacal dust model can account for these 
excesses. We emphasize the approximate nature of this 
estimate and underscore the point that the actual dust 
configuration may be far more complicated. Neverthe- 
less, a single-ring approximation yields a radial location 
as indicated in Fig. the probability distribution for 
the inner ring radius. The uncertainty encloses the 66% 
confidence interval in the probability distribution. 

It is evident in Fig. [3] that the simple power-law de- 
scription for cr(r) in the outer ring fails beyond the outer- 
ring edge. We expect this discrepancy to have no effect 
on our estimate of the zodiacal dust component. How- 
ever, we explore it below with a more detailed model to 
better understand outer-ring structure. 

3.2. An Outer Ring with Two Components 

To date, the wide variety of published imaging and 
photometric measurements of HR 4796A have not been 
used simultaneously to constrain model simulations. An 
oft-used single-ring model fits some of the individual im- 
ages plausibly well. Here, however, we show that a more 
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complicated model of the outer ring is required to fit all 
available data. We first find the single-ring model that 
provides the best simultaneous fit to all data, includ- 
ing Keck/MIR images, the HST-NICMOS/N1R image, 
and all extant flux density measurements of HR 4796A 
(Table, f). We assume flux values for the zodiacal con- 
tribution that were estimated in the previous section; 
since these contribute negligibly to the other data, only 
outer-ring parameters are varied. For the coronagraphic 
HST/NICMOS image, a circular region of radius 0.6" 
centered on the star is masked, so that only the ansae 
of the outer ring are considered. This accords with S99 
who point out that emission inside this region is largely 
an artifact of photosphcric light scattered off the inner 
edge of the coronagraphic hole. In image fitting, models 
are freely scaled and thus only fit to the intensity pattern 
of the data. The total fluxes at the imaging wavelengths 
are fit as part of the SED. 

To minimize the number of free parameters, we make 
one additional simplification at the onset of detailed 
modeling. Searches for best-fit images revealed that 
the NE ansae should have an optical depth that is 20% 
greater than the SW ansae's optical depth. So we incor- 
porate this asymmetry into all our models by introducing 
a sinusoidal azimuthal variation into the optical depth 
description, such that the function has a crest at the NE 
peak and a trough at the SW peak. The amplitude of 
this azimuthal variation is kept constant across all model 
fits, thus ridding us of one free parameter. 

The inadequacy of the single-ring model is most ev- 
ident in comparison with images at A = 12.5 /im and 
1.1 /im and in a discrepancy between model and observed 
fluxes at thermal infrared wavelengths. In Fig. Etc), the 
best-fit single-ring model fails to simulate the extended 
emission in the 12.5 /im image (Fig.JSJa))- It is also dis- 
crepant with the MIR flux densities as evident in Fig. 
Finally, emission profiles at 1.1/zm (solid line) in Fig. 
are not well matched by a single-ring model (dot-dashed 
line) fit to all data. The peaks of model emission are 
offset toward the star, and model flux is too high at the 
outer radius and too low at the inner edge. This mis- 
match is consistent with the positional offset between 
peaks of scattered- light and thermal emission in Fig. 6(b) 
of TOO, an overlay of the HST/NICMOS countour map 
on a 18.2/xm Keck/OSCIR image of HR 4796A. However, 
it was unclear if this effect was simply due to the higher 
temperature at the inner edge of a single outer ring. Ev- 
idence outlined above indicates the discrepancy in peak 
locations is also a consequence of a more complex dust 
distribution. 

Some hints as to how the model can be improved are 
evident in fits to some of the images. A single uniform- 
density (7 = 0) ring model fit to the MIR images alone 
gives outer-ring widths of ~ 60 AU and grain radius ~ 
5 /im. This is close to the grain radius estimates of TOO, 
where the chief constraint is the relative brightnesses of 
the outer disks at 10 and 18 /Ltms. On the other hand, fits 
to the SED alone yield a 180-AU-wide ring with roughly 
16 /im grains. Smaller disk sizes, yield larger grain sizes. 
In fact, the Wien side of the outer disk's SED demands 
small grain sizes, while the Rayleigh- Jeans side of the 
outer disk's SED (the submm flux densities) demand 
large grain sizes. Independant fits to the HST/NICMOS 
image yield a narrower ring, however, ~ 17 AU. These 



TABLE 2 

The table below shows the unreduced x 2 from fits to 
each data set for the single-outer-ring model and the 
2-component-outer-ring model. the top row of the table 
shows the number of data points in each data set. the 
next two rows show the corresponding x 2 s for each 
model. the single-outer-ring model (1-comp) had 9 
parameters, while the 2-comp had 13 parameters. 
According to the Bayesian Information Criterion, which 
takes into account the total unreduced \ 2 for- each 
model and the number of data points and free parameters, 
the 2-comp model is statistically preffered to a 
significant degree (see details in text). the x 2 values 
have been rounded to the nearest integer. 



Model 24.5/im 20.8/tm 12.5/im 1.1/an SED total 

No. data pts. 13 12 12 53 16 106 

2-comp x 2 13 15 14 53 15 110 

1-comp x 2 19 16 22 70 27 154 



discrepancies suggest that the outer ring has more than 
one component, and that it might be better approxi- 
mated as two rings, each with a different characteristic 
grain size. 

Comparisons of the 2-component outer-ring model 
with the single-outer-ring model and data are shown in 
Figs. El & □ The entire SED from the best-fit 2- 
component model is shown in Fig|SJ In all cases sig- 
nificant improvements are evident. However, for a more 
robust quantitative analysis of which model is preferred, 
we can employ the Bayesian Information Criteria [BIC) 
(Liddle 2004; Mukherjee et al. 1998; Jeffreys 1961). This 
is a model selection criteria which penalizes models for 
having extra parameters: BIC — —2lnL + klnN , where 
L is the likelihood function e~ 1 / 2S " x « , k is the number of 
free parameters and N is the number of data points. The 
model with the lower BIC is preferred. A difference of 
2 between BICs indicates positive evidence against the 
high-BIC model, while a difference of 6 indicates strong 
evidence. The top row in Table. 2 shows the number of 
data points for each of the data sets (images and flux 
measuremtents) . The next two rows show the unreduced 
X 2 s (rounded to the nearest integer) for each of the data 
sets from the two best-fit models, k for the 1-comp model 
is 9, while k for the 2-comp model is 13. From these num- 
bers, the BIC for the 1-component model is 196 and the 
BIC for the 2-component model is 171. The magni- 
tude of this difference provides strong evidence that a 
2-component model is mandated by the data. 

The best-fit model (2-component outer-ring) is shown 
at MIR wavelengths in Fig. |5J Single-ring best-fit pa- 
rameters are: R in = 64 ± AAU , AR = 80 ± 6AU, 
7 = 4.88 ± 0.2, A = 31 ± 3 /im (a = 21 /im), a = 
3.3 ± 0.2 x 10" 3 and albedo= 0.3 ± 0.02. Uncertain- 
ties were estimated through the Bayesian approach. The 
2-component outer-ring model parameters are: R% n \ — 
45 ± 5AU, Ai?i = 80 ± 15AU, X 01 = 5 ± 2 /im (a 2 = 3.3 
/im), a i = 2.1 ± 0.2 x 10~ 3 (implying a total cross- 
sectional area of 90 AU 2 ) and albedo, u>i — 0.35 ± 0.03, 
R m2 = 66 ± IAU, AR 2 = 14 ± IAU, A 02 = 38 ± 5 
/im (a = 25 /im), a 2 = 3.8 ± 0.25 x 10~ 2 (implying 
a total cross-sectional area of 245 AU 2 ) and albedo, 
u>2 = 0.18 ± 0.02. The 7 parameter was not used for 
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Fig. 5. — Comparison of 12.5 fim image (a) with simulations from 2-component outer-ring model (b) and a single-outer-ring model (c). 
Both simulations are from best-fit models to all extant data. Contour levels are same as in Fig. UTal 
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Fig. 6. — Strip cuts through the l.lfim HST/NICMOS image are shown as solid lines. Cuts through a simulated image from the 
2-component outer-ring model are shown as dashed lines. The dash-dot-line represents the single-component model. All 8-AU-wide strip 
cuts are along a PA of 26° and pass through both NE & SW peaks of emission. The NE portion is shown on the left panel and the SW 
portion is shown on the right panel. The section in between was masked by the NICMOS coronograph. 



either ring. The probability distributions, obtained by 
the Bayesian approach, for some of these parameters are 
shown in Fig. 1101 

4. DISCUSSION 



Analysis of multi-wavelength HR 4796A images and 
flux densities yields a refined estimate of the properties of 
its circumstellar dust as follows. Radial structure of the 
outer ring is best approximated by two components: a 
narrow ring of ~50 /im grains between 66 and 80 AU and 
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Fig. 7.— Flux densities of HR 4796A at A = 18.2, 20.8, and 
24.5/iin from Table. 1, shown with error bars. The dash-dot line 
represents the total SED from the 2-component outer-ring model. 
The dotted line represents the total SED from the single outer-ring 
model. 



a surrounding wider ring of ^7 /im grains stretching from 
45 AU to about 125 AU. The presence of excess emission 
at the stellar position is confirmed. It is unresolved by a 
0.37" FWHM PSF (corresponding to a 12.5 AU radius at 
the 67 pc distance of HR 4796A) and has a temperature 
of ~ 260 K. Bayesian parameter estimation suggests this 
dust is ~ 4 AU from the star with a total grain cross- 
sectional area of 0.055 AU 2 . 

4.1. Exozodiacal Dust 

Evidence for warm exozodiacal dust was first identi- 
fied for HD 98800 (Zuckerman & Bccklin 1993). This 
K-dwarf quadruple system of two spectroscopic binaries 
(Torres et al. 1995) exhibits one of the strongest dust 
signatures of any in the IRAS catalog. Its infrared ex- 
cess peaks at 25 /xm and is associated entirely with the 
northern spectroscopic binary, HD 98800 B (Koerner et 
al. 2000). Modeling of the spectral energy distribution 
of the dust suggests a ring of particles in the few-micron 
size range with an inner radius of a few AU and an outer 
radius that may stretch will into a giant-planet zone. 
Subsequent detections of warm excess from the main- 
sequence A3 star, £ Leporis, have also been intepreted 
as orginating from collision in an asteroid belt (Chen & 
Jura 2001). 

The origin of a warm component to circumstellar dust 
around HR 4796A may be similar to that of our own Solar 
System. Zodiacal dust is supplied predominantly by the 
asteroid belt, Jupiter Family comets (JFCs), long-period 
comets, and perhaps by dust from the Kuiper-belt and 



the Oort Cloud (Flynn 1994, Dermott et al. 1992, Durda 
& Dermott 1992, Liou & Dermott 1993). At 1 AU, about 
45% of the dust in the ecliptic plane is believed to come 
from the asteroid belt and JFCs. Within a sphere of 
radius 1 AU, however, 89% of the dust is from comets. 
The fractional surface density of zodiacal dust grains is 



<j(r) 



-1.45 



(Hahn et al. 2002) with grain sizes as 



large as 100 /im (Grogan et al. 2001). 

An asteroid belt may be the source of exozodiacal dust 
grains around HR 4796A. Their total grain cross-section 
of ~ 0.055 AU 2 and effective grain size of 28 /im yields 
a minimum mass of ~ 10 -3 Mq, under the assumption 
that the grain density is given by p=2500kg/m 3 (Jura 
1998). Ice sublimation occurs instantaneously at 4 AU 
from the star (Isobe 1970), so these grains are almost 
certainly refractory. The blow-out size due to radiation 
pressure is 7 /im (calculated according to Backman & 
Paresce 1993, Gustafson 1994) and indicates this is not 
an effective dispersal mechanism for the grains we de- 
tect. Due to relatively high fractional surface density 
(<r = 0.0046), the collisional lifetime is shorter than the 
Poynting-Robcrtson lifetime even at 4 AU. Grains are 
primarily being destroyed by collisions, but probably also 
originate in a cascade of collisions from larger bodies. As- 
suming a grain size distribution, n(a) ~ a -3 5 , and 1000 
km planetesimals for the largest bodies, the total mass 
of orbiting material is ~ 0.6 Mq. 

Another possible source for HR 4796A's exozodiacal 
dust is a more distant asteroid or exo-Kuiper belt. Grains 
may spiral in toward the star under the influence of 
Poynting- Robertson drag but be trapped by resonant in- 
teractions with terrestrial planets. The dust may come 
only from more tenuous parts of the disk where collisional 
destruction is not the dominant dispersal phenomenon. 
For our Solar System, however, Liou et al. (1996) showed 
that Kuiper belt grains in the range 9 to 50 /zm are un- 
likely to survive collisions with interstellar grains (as op- 
posed to interplanetary grains) on their way to the Sun, 
and 80% of all migrating grains are likely to be thrown 
out of the Solar System by interactions with the giant 
planets. However dust can be trapped in planets inte- 
rior to an asteroid belt, and in these cases grains roughly 
30 /im or greater in size seem to be the best candidates 
(Jackson & Zook 1992). Large asteroid grains can even 
be injected into interior resonances with exterior planets 
due to the initial outward thrust of radiation pressure 
upon release from a parent body. As they leave these 
resonances they usually have lower eccentricities and or- 
bital inclinations that render them susceptible to capture 
in resonant orbits with interior planets. 

Comets that travel close to the star may provide a 
source of exozodiacal dust grains around HR 4796A. The 
comctary component of zodiacal dust is thought to be 
highly variable according to simulations done by Napier 
2000 and may increase the dust creation rate from a 100 
to 1000 fold within periods of 100,000 to 10,000 yrs. Par- 
ticles that are most likely to be trapped in our own Solar 
System have sizes between 10 and 100 /im (Jackson & 
Zook 1992). Dust with low eccentricities and low orbital 
inclinations are again more likely to be captured in reso- 
nant orbits. It should be noted that while the radiation 
pressure blow-out size for the Solar System is roughly 0.9 
/im, it is 7 /im for HR 4796A. The optical depth of the 
zodiacal dust at 25 /im is roughly 10~ 7 when measured 
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Fig. 8. — Flux densities of HR 4796A from Table 1. Wide-band photometry ignored where narrow-band photometry available. The 
dashed line represents the total SED from the best-fit 2-component outer-ring model. The solid line is the contribution from the photosphere 
(derived from Kurucz 1993 model). The dotted line is the contribution from the wide outer-ring component. The dash-triple-dot line is 
the contribution from the narrow outer-ring component. The dash-dot line is the contribution from the inner zodiacal dust component. 
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Fig. 9. — The simulated images from the best- fit 2-component outer-ring model at the same wavelengths as in Fig. Q (a) Simulation 
of emission at 12.5 /im. (b) Simulation of emission at 20.8 fim. (c) Simulation of emission at 24.5 fim. The contour levels for each figure 
are the same as those in Fig. HI 



along the ecliptic (Spiesman et al. 1995). For HR 4796A, 
assuming similar dust dispersal times scales, this would 
imply a dust production rate roughly 10 4 times higher 
than expected for the Solar Sytem. Grain velocities also 
differ in the two systems. Thus it is best to undertake 



system-specific simulations before speculating too much 
about the origin of exozodiacal dust grains. Regardless 
of origin, the main explanation for the persistence of dust 
in the exozodiacal region is probably resonance trapping 
by interior planets as found in the Solar System. 
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Fig. 10. — Probability distributions for selected 2-componcnt outer-ring model parameters. Most probable values correspond to peaks 
in the probability distributions. The uncertainties are estimated as the 66% confidence intervals (where our distributions arc binned 
finely enough) or the shortest range of parameter values that encompasses 66% of the total probability. The most probable values and 
uncertainties for these parameters are detailed in the text. 



4.2. Properties of HR 4796 A 's Outer Ring 

A cold, massive outer ring of dust particles is more 
readily detected around stars other than the Sun. Ex- 
amples other than HR 4796A include Fomalhaut (Hol- 
land et al. 1998); e Eri (Greaves et al. 1998), and Vega 
(Holland et al. 1998; Koerner et al. 2001; Wilner et al. 
2002). However, such a structure is more difficult to de- 
tect in our own Solar System where confusion from the 
inner zodiacal dust interferes with emission from outer 
dust grains (Backman et al. 1998). The discovery of a 
sizeable aggregation of Kuiper Belt Objects (see review 
by Luu & Jcwitt 2002) implies an associated population 
of dust grains (Backman et al. 1995). Larger grains are 
likely to be better retained in the Kuiper Belt region 
by mean motion resonances while smaller grains diffuse 



under the influences of P-R drag, solar-wind drag, and 
radiation pressure (Holmes et al. 2003). 

The HST/NICMOS image of HR 4796 A at 1.1 /zm es- 
tablishes that most of the scattered light around HR 
4796 A arises within a 17 AU ring centered at 70 AU 
(S99). A wider low-density component is also needed to 
fit Keck/MIRLIN images of thermal emission. The spec- 
tral energy distribution can be fit by a single narrow ring, 
but only if the large-grain contribution to the size distri- 
bution is boosted significantly. LL03 used n(a) ~ a~ 2,9 
to do this, but their models fail to match the surface 
brightness distribution in thermal infrared images. In- 
clusion of two distinct effective grain sizes resolves an 
apparent inconsistency in the interpretation of different 
data sets. Thermal infrared images modeled by TOO re- 
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quired small grains (~ 3 /irn), but a fit by K98 that 
included an upper limit to the flux density at 800 /im 
required the presence of large grains (~ 30 fim) to fit 
the Ray leigh- Jeans side of the SED. The latter result 
is confirmed here by fits which include more recent de- 
tections of emission at 450 and 800 (im (Holland et al. 
1998) and a new measurement at 350 fim. A two-part 
ring with separate grain-size distributions is necessary to 
reproduce observations at all wavelengths. 

Properties derived here for HR 4796A's outer ring are 
consistent with theoretical predictions for both our Solar 
System's outer dust and for exo-Kuiper grains in gen- 
eral. The collision of large planetesimals is expected to 
produce a cascade of smaller-sized particles with number 
distribution dn(a)/da ~ a -3 ' 5 (Dohnanyi 1969). Sub- 
sequent grain evolution is affected largely by radiation 
pressure; Poynting-Robertson drag is a relatively small 
effect. Radiation pressure effects can be parametrized by 
0: 

(3(D) = (115O/pD)(L*/L )(M /M*) 

where D is grain diameter in /ims, and p is grain density 
in SI units (Gustafson 1994, W99). Appropriate values 
for HR 4796A are 2.5 M (Jura 1998) and p=2500 

kg/m 3 and yield (3(D) = 3.5/ D. Large particles with 
D > 35 correspond to (3 < 0.1 and will remain in original 
orbits close to parent bodies. Grains in the size range 35 
/j,m > D > 7 fim are (3 critical with values 0.1 < (3 < 0.5. 
These experience orbital evolution which extends their 
distribution both inwards and outwards. The smallest 
grains (D < 7 /im and (3 > 0.5) have been called "(3 me- 
teoroids" and are blown out of the system on hyperbolic 
orbits. 

This picture of grain evolution provides a natural ex- 
planation for the properties of the outer ring around HR 
4796A. In this scenario, a large-grain ring centers on a 
system of planetesimals surrounded by an outer region 
of smaller grains diffused under the influence of radiation 
pressure. Large grains in the narrow ring have (3 < 0.1, 
D > 35 /Urn or D ~ 50 ^m in our model, so they stay close 
to parent bodies. Smaller grains are 11 f3 critical" or "(3 
mctcoroids," since their size is so close to 7 pm (D ~ 6.6 
jitm) threshold. However, "(3 critical" is more likely, since 
u (3 meteoroids" have a "blow-out" lifetime of order 100 
yrs, while "(3 critical" dust has a collisional lifetime of or- 
der 1000 yrs. This idea is supported by comparison with 
expected size distributions. The large- and small-grain 
populations have a diameter ratio of about 7:1 and cross- 
sectional areas 90 AU 2 for small grains and 245 AU 2 for 
large grains. So the number of grains is approximately 
in the ratio 20:1 for small:large grains. In contrast, the 
initial ratio for a collisionally generated size distribution 
should be (7/50)" 3 5 - 1000 : 1. This underabundance 
of small grains agrees well with the hypothesis that many 
"/? meteoroids" have already left the system. 



4.3. Furthur Work 

The modeling approach presented in this work was 
designed to determine the simplest disk morphology re- 
quired by observations with a minimum number of under- 
lying assumptions. The resulting dust distribution sug- 
gests significant physical interpretations. These should 
be tested with a more complex analysis that uses a more 
realistic representation of dust grain properties, includ- 
ing physically meaningful size distributions and plausible 
assumptions about grain composition. The latter should, 
in turn, be coupled to models that simulate the dynam- 
ical evolution of dust grains under the influence of plan- 
etary bodies in the HR 4796A system. 

Further progress in observing HR 4796A's radial 
dust distribution will be difficult until a next genera- 
tion of large- aperture telescopes comes on line, since 
Keck/MIRLIN images are already close to the diffrac- 
tion limit of a 10-m telescope. ALMA may provide sub- 
millimeter images with resolution improved by a factor 
of 5. The next generation telescopes like the 30 meter 
"Giant Segmented Mirror Telescope" , the "Thirty Meter 
Telescope" and the 100 meter "Overwhelmingly Large 
Telescope" will have also improved resolution and greatly 
enhanced sensitivity, and high-dynamic range techniques 
that are designed to search for planets may provide more 
detail on reflected light. Meanwhile, progress in observa- 
tions of exozodiacal and exo-Kuiper dust configurations, 
may be better advanced by studying systems that are 
closer than the 67 pc distance to HR 4796A. The en- 
hanced sensitivity of Spitzer Space Telescope will iden- 
tify nearby debris-disk systems for which spatial reso- 
lution will be improved. Nulling interferometers (Keck 
and LBT) will also bring enhanced capability to study 
their inner dust regions. Ultimately, coordinated planet 
searches will be able to verify the relationship between 
circumstellar dust signatures and the presence of a plan- 
etary system. 

Portions of this work were carried out at the Jet 
Propulsion Laboratory, operated by the California Insti- 
tute of Technology under a contract with NASA. Data 
presented herein were obtained at the W.M. Keck Ob- 
servatory, which is operated as a scientific partnership 
among the California Institute of Technology, the Uni- 
versity of California and the National Aeronautics and 
Space Administration. The Observatory was made pos- 
sible by the generous financial support of the W.M. Keck 
Foundation. The authors wish also to recognize and ac- 
knowledge the very significant cultural role and reverence 
that the summit of Mauna Kea has always had within the 
indigenous Hawaiian community. We are most fortunate 
to have the opportunity to conduct observations from 
this mountain. 
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